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Abstract. We propose to apply atom-chip techniques to the trapping of a single atom in a circular Rydberg 
state. The small size of microfabricated structures will allow for trap geometries with microwave cut-off 
frequencies high enough to inhibit the spontaneous emission of the Rydberg atom, paving the way to 
complete control of both external and internal degrees of freedom over very long times. Trapping is achieved 
using carefully designed electric fields, created by a simple pattern of electrodes. We show that it is possible 
to excite, and then trap, one and only one Rydberg atom from a cloud of ground state atoms confined 
on a magnetic atom chip, itself integrated with the Rydberg trap. Distinct internal states of the atom are 
simultaneously trapped, providing us with a two-level system extremely attractive for atom-surface and 
atom-atom interaction studies. We describe a method for reducing by three orders of magnitude dephasing 
due to Stark shifts, induced by the trapping field, of the internal transition frequency. This allows for, 
in combination with spin-echo techniques, maintenance of an internal coherence over times in the second 
range. This method operates via a controlled light shift rendering the two internal states' Stark shifts 
almost identical. We thoroughly identify and account for sources of imperfection in order to verify at each 
step the realism of our proposal. 
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^ 1 Introduction 



In recent years we have witnessed convergence between 
many fields of physics previously considered disparate. 
This has been particularly true in the case of solid state 
and atomic physics. On one side, the possibility of design- 
ing nanometric devices and the understanding of quantum 
phenomena in devices such as Josephson junctions has al- 
lowed experimentalists to control the quantum coherence 
of semi- pQ or super-conducting [2113] systems, exactly as 
was done many years ago for atoms or molecules. On the 
other, atomic physicists have been able to produce Bose- 
Einstein condensates 0] and degenerate Fermi gases in 
quasi ideal situations, wonderful tools for probing many 
body theory and mesoscopic physics. 

Atom-chip experiments jHHZj open the way to inves- 
tigations at the frontier of these two fields. Aside from 
their relative ease of use and potential as atom inter- 
ferometers, they provide an ideal environment in which 
atomic ensembles can be integrated with, and coupled to, 
devices on a microcircuit. Moreover, they allow for the 
study of atom-surface interactions with a completely new 
degree of control. Separately, experiments on circular Ry- 
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dberg atoms have proven them extremely sensitive probes 
of electric and magnetic fields, both static and dynamic, 
making ideal tools for the investigation of numerous quan- 
tum phenomena, atom-surface interactions included. This 
remarkable sensitivity, in conjunction with their relative 
ease of detection, makes a single circular Rydberg atom 
an excellent probe of microwave field intensities ranging 
from one to a few tens of photons |5]. The coupling of a 
single such atom to a high Q superconducting microwave 
cavity has allowed our group to study, for example, the de- 
coherence of mesoscopic ensembles of photons [0] and to 
realize elementary quantum logic operations In these 
experiments, circular Rydberg atoms are excited from a 
thermal atomic beam (mean velocity 350 m/s) crossing 
a cryogenic setup in which the superconducting cavity is 
mounted. The use of thermal atoms intrinsically limits 
the maximum interaction time between the atom and the 
microwave field. Moreover, in order to reduce to a neg- 
ligible level the proportion of unwanted events in which 
we have more than one circular Rydberg atom present at 
once in the cavity, we are forced to work with a very low 
rate of excitation towards the Rydberg levels. Aside from 
increasing dramatically the acquisition time, this lack of 
a deterministic atom source is an inconvenience from the 
point of view of quantum information processing. 
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The integrated trapping and manipulation of circular 
Rydberg atoms on a chip proposed here would lift these 
limitations. If we are to take advantage, however, of the 
longer interaction times achievable with such trapped Ry- 
dberg atoms, we must also combine in the same setup the 
capability for the extension of the natural lifetime and long 
term coherent manipulation. For short term manipulation 
over times of the order of the millisecond, this being much 
shorter than their lifetime in free space, one could simply 
excite the Rydberg atoms from an initial sample of cold 
trapped atoms and study them in free fall JT] . In Ref. \F2\ 
we proposed an electric trap able not only to store an indi- 
vidual circular Rydberg atom, but also to inhibit its spon- 
taneous emission and maintain an internal coherence. In 
this paper we present in much more detail how we assess 
the trap's performance, how to integrate it on a chip and 
how we simulate the different imperfections of the setup. 
We not only present the trap geometry of Ref. ^2 but 
also show that it is possible to extend the proposal to a 
smaller design where the atom is brought closer to the chip 
surface, a design therefore better adapted to atom-surface 
interaction studies. The principle of the trap, along with 
its implementations in realistic geometries, is discussed in 
Sec. We show in Sec. 01 that the excitation of a sam- 
ple of ground state atoms, initially trapped on a magnetic 
atom chip, can lead to the preparation of a single Rydberg 
atom, taking advantage of the dipole-blockade effect |13| . 
As explained above, the maximum duration of an experi- 
ment with this atom is not limited by the storage time but 
by the lifetime of the circular Rydberg state. We discuss in 
Sec.rjjthe efficiency of the spontaneous emission inhibition 
when we place our atom chip within a structure exclud- 
ing the millimetre-wave decay channel. We show that the 
radiative lifetime could realistically be pushed into the sec- 
ond range. We present finally in Sec. 0a technique taking 
advantage of a controlled light shift of the Rydberg levels 
allowing for the maintenance and control over a similar 
duration of the coherence of an atom in a superposition 
of two trapped states. 

2 A microfabricated electrodynamic trap for 
Rydberg atoms 

The large coupling of Rydberg atoms to static or time- 
varying fields presents us with numerous possible trapping 
techniques. It has been proposed, for example, to use pow- 
erful, far-detuned laser beams to create very tight traps 
for Rydberg states. These traps exploit the ponderomotive 
force [ID experienced by the atom due to the fast oscil- 
lation of its valence electron in the laser field. It remains 
hard however to fulfill all the criteria allowing for coherent 
manipulation of Rydberg levels in any such trap. For in- 
stance, one can see that the very large laser power required 
by the ponderomotive trapping is hardly compatible with 
the cryogenic environment necessary to avoid blackbody 
radiation in the microwave domain |8j and the consequent 
destabilization of the Rydberg levels. 

Another possible solution is to use d.c. magnetic fields 
and conventional atom-chip architectures. The magnetic 



moment of a Rydberg level (principal quantum number n) 
can be of the order of n/is, where hb is the Bohr mag- 
neton. The same atom-chip wires could be used for the 
trapping of ground and Rydberg state atoms. One possi- 
ble problem is that the Rydberg state might be affected by 
its interaction with still-trapped ground-state atoms in its 
vicinity. Moreover, the transition frequency between dif- 
ferent Rydberg states is strongly perturbed by the Zeeman 
effect. We could not find a way to suppress the resulting 
random dephasing, in contrast to that due to the Stark 
effect, which can be partially compensated as explained 
in Sec. El 

In the following we therefore use time-varying electric 
fields to trap the Rydberg atoms. 



2.1 Circular Rydberg atoms in an electric field 

We consider the case of Rydberg atoms in circular states. 
These have a large principal quantum number n, of the 
order of 50 in our discussion, along with maximal angu- 
lar and magnetic quantum numbers (I = \m\ = n — 1). 
The atom therefore has both a large energy and a large 
angular momentum. The classical analog of the electron's 
wavefunction is a Keplerian circular orbit around the nu- 
cleus. To be concrete we will consider 87 Rb atoms. To a 
very good approximation one can neglect the fine struc- 
ture corrections to the circular levels' energies and assume 
that they are equal to those of Hydrogen. We will label \g) 
the circular Rydberg state with principal quantum num- 
ber n = 50 and |e) that with n = 51. In the presence of an 
electric field the good quantum numbers are n, m and the 
parabolic number m [El. The circular state corresponds 
to the case \m\ = n — 1, ni = Q (it is simultaneously an 
eigenfunction in the {n, l,m} and {n,ni,m} representa- 
tions) and the normal to the plane of its circular orbit 
aligns parallel to the electric held. Due to the Stark effect 
the degeneracy of the manifold of levels of equal n but 
differing n\ and m is lifted. Under these conditions the ra- 
diative lifetime of the levels \g) and |e) is extremely long, 
r sp ss 30 ms. For a given electric field amplitude E the 
energy of each level can be calculated using either a direct 
diagonalization of the Stark Hamiltonian (see Sec. I5.r?fl or 
a perturbative approach. Using the latter method, to sec- 
ond order, the energy £{E) of a given level \n, m,m) is 
given by: 

£(£)=£'">+ £<'>+£' 2 >, (1) 

£ "" = -5?' (2) 

£ (1) = -^(n- 2ni - \m\ - l)nE, (3) 

£ {2) = - — [17n 2 -9m 2 + 19 
16 

-3(n- 2ni - \m\ - l) 2 ]n 4 £ 2 , (4) 

where all quantities are expressed in atomic units (Energy: 
£ a u = 4.360 x 10~ 18 J, electric field: E a . u . = 514.2 x 
10 9 V/m). In the presence of a time- varying electric field, 
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Fig. 1. Level spectrum in presence of an electric field for the 
multiplicities n — 50, 51 and 52. Only the levels with the largest 
angular momentum are presented (m > n — 3). A 7r-polarized 
microwave field will couple levels of equal m, for example those 
of m = 49 separated from the others by a box in this figure. 
A (j-polarized microwave field will, in contrast, couple levels of 
m differing by ±1, for example levels \g) and \e). 



the normal to the plane of the orbit follows the field direc- 
tion so long as the characteristic frequency of the change 
in electric field remains very small compared to the tran- 
sition frequency ^ stark = £^/h between the neighbouring 
non-circular levels. In an electric field of 400 V/m, typical 
of the fields considered in this article, this frequency is 
of the order of 400 MHz, much higher than all the trap- 
ping frequencies found later, and therefore than the rate 
of change in angle of the field. We can therefore in all the 
following consider the evolution in the time-varying field 
to be adiabatic and that the atom remains in the circular 
Rydberg state [TS] . 

Fig. n presents the energy levels in the presence of an 
electric field for the principal quantum numbers n rele- 
vant to our discussion. It is important to note that cir- 
cular states, having no permanent mean electric dipole, 
experience no linear Stark effect (£^ — 0). They are, 
however, much more highly polarizable than ground state 
atoms and have an accordingly large quadratic Stark ef- 
fect [£i » (XiE 2 for i = (g,e), where Q!j is half the polar- 
izability]. Its value, a g = -203,2 Hz/(V/m) 2 for n = 50 
and a e = -228,7 Hz/(V/m) 2 for n = 51, is 9 orders of 
magnitude larger than that of the ground state of Hydro- 
gen. Nonetheless, circular states are high-field seekers and 
cannot therefore be trapped by any configuration of d.c. 



electric fields, a maximum of the electric field modulus in 
vacuum being forbidden by Maxwell's equations. 



2.2 Trapping 

A similar situation is encountered in the case of charged 
particles, such as ions. In that particular example, Max- 
well's equations prevent us from obtaining an extremum 
of the electric potential in a region devoid of charge. It is 
nevertheless possible to trap ions if one uses a.c. electric 
potentials, as in the case of Paul traps JTj. Building on 
these ideas it has been shown by Peik ^B] that polarizable 
atoms can be trapped if one combines a relatively strong, 
d.c. and homogeneous electric field E\ — E\U Z (where u z 
is a unit vector along the vertical) and an inhomogeneous 
a.c. field E%(r,t) deriving from a time- varying hexapo- 
lar potential. As a first step we consider the case where 
£3 <C Ei. The field is therefore almost homogeneous over 
the trapping region and small departures from E\ are re- 
sponsible for the confinement. Let V\{r) and V%(r,t) be 
the potentials associated to E\ and E 3 (r,t) respectively. 
One can write these potentials in the form: 



Vi(x,y,z) = 
V 3 (x,y,z,t) = 



Uiz 
2z 3 



(2z 3 - 3zx 2 - 3zy 2 ), 



(5) 
(6) 



where the Ui are quantities homogeneous to electric po- 
tentials and zq is a length of the order of the electrode size. 
If one considers a circular Rydberg atom following adia- 
batically the electric field variations, its potential energy 
in the quadratic Stark approximation is given by: 



£ = a\E x - 



2E 1 



Ei 



(7) 



Assuming E% <C E\ and neglecting the constant terms in 
the potential energy we arrive at: 



£ 



2aE l E 3 ■ u z 

U x dV 3 {x,y,z,t) 



2a 



z dz 



v 2 )- 



(8) 



We see that we now have an expression for the potential 
energy quadratic in x, y and z, as in conventional ion 
traps. It is important to note that if we have trapping 
along directions x and y (i.e. UiU^t) > for a < 0) then 
we necessarily have antitrapping along direction z, and 
vice versa. Moreover, if U3 varies sinusoidally [i.e. t/3 (t) = 
U30 cos(wf)] then the equation of motion along Ox can be 
written: 



d 2 x QaUiUw 
mRb ^u2 4 cos(Lut)x = 0. (9) 



Equation|^lcan be transformed into the well known Math- 
ieu form: 
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fx 
dr 2 



2q x cos(2r)a; = 0, 



(10) 



where r = wt/2 et q x = ^aUiUzo/mRbUj 2 z$. 

The equation of motion is of the same form along the 
other directions with q y = q x and q z = —2q x . There exists 
a stable solution, confined in space, if < 0.907, i = 
x,y,z. 

In order to compensate the force of gravity, taken to 
be along —z, it is also possible to add to V\ and V3 a d.c. 
quadrupolar potential V 2 (r) = C/ 2 (x 2 + y 2 — 2z 2 )/zq. If 
U2 <C Ui, the same analysis as before, considering only 
the terms in U\ and U2, leads to: 



2a 



Ui dV 2 (x,y,z) 



dz 



-4a — o — z, 



compensating gravity if U2 = ragz^/ AaU\. 



(11) 



2.3 Proposed geometries, calculation of potentials and 
fields 
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It is of course possible to design a geometry of electrodes 
which perfectly maps the desired hexapolar geometry [see 
Fig. 12(a)]. It requires two rings and two end caps. How- 
ever, being closed in all three dimensions, such a geometry 
is hardly compatible with efficient atom loading, nor with 
the idea of building all the trapping components on a mi- 
crofabricated structure. We propose instead a design (trap 
A) based on two chips facing each other [see Fig. 0(b)]. 
Each chip could be made by standard lithographic meth- 
ods and consists of two electrodes: a disk, playing the role 
of end cap, surrounded by a plane. The diameter of the 
disks is 1 mm and the outer electrodes would extend to 
y — ±1 cm and even further in x, the electrodes forming 
part of a long waveguide (see Sec. I5.2fl . To ease the nu- 
merical calculation of the potentials we cut the electrodes 
off at p = \J x 2 + y 2 = 2 mm and confirmed afterwards 
that this had a negligible effect upon the potential within 
the trapping region. 

The voltage applied to each electrode is the sum of 3 
terms: 

— A static voltage ±Ui creating a homogeneous field E\ 
along Oz, as in a plane capacitor. 

— A static voltage U2 on the outer ring of both chips cre- 
ating a potential with a large quadrupolar component 
of the same form as V^ir). This component, in con- 
junction with Ei as seen above, will create a force on 
the atom compensating gravity. 

— A time- varying voltage U3 (t) creating a potential with 
a large hexapolar component of the same form as V3 (r , t) 
This component, in conjunction with E\ as seen above, 
allows for trapping. 

It is, however, important to note that the electric po- 
tential associated to U2 (resp. E/3) is clearly not an ideal 
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Fig. 2. For all figures there is a cylindrical symmetry around 
Oz and the electrodes are shaded according to the phase of 
the oscillating potential 173 . (a) Electrode geometry creating 
an exact hexapolar potential, (b) Section of trap A in a ver- 
tical xOz plane, with applied potentials. The diameter of the 
inner electrode is 1 mm. The plate spacing, 1 mm, is appro- 
priate for spontaneous emission inhibition. The analogy with 
the geometry of (a) is conspicuous, (c) Geometry of trap B. 
The lower plane is not drawn to scale, the inner disk electrode 
having a diameter of 100 /jm and the ring around it being of 
width 100 fj,m. The plate spacing remains 1 mm. 



quadrupole (resp. hexapole). The resulting electric fields 
are therefore more complicated. As an illustration, if the 
factor 77 of Fig. (b) is set equal to 1, U3 will create a 
non-zero electric field along Oz at the centre of the trap 
O. This will create a time- varying electric field, and there- 
fore atomic energy, at trap centre, resulting in heating. A 
correct choice of 77 allows us to better approximate a per- 
fectly hexapolar potential, and cancel this electric field. 



John Mozley et al.: Trapping and coherent manipulation of a Rydberg atom on a microfabricated device: a proposal 5 



More specifically, our final goal is to be able to calcu- 
late the electric field at any position inside the trap and 
at any time. Moreover, it is important to assess edge ef- 
fects as well as the effects of the finite gap between the 
electrodes. We have therefore opted for a numerical ap- 
proach. We numerically calculate the potential inside the 
trap on a grid of points using the software SIMION |19) . 
In order to limit numerical errors in the calculation of the 
gradient we then make a fit of the potential on the basis 



of the first seven spherical harmonics Y{ 



centred 



on O. Due to the cylindrical symmetry around Oz only 
the m = terms are present. Thanks to the linearity of 
Maxwell's equations one can add separately the contribu- 
tion of each voltage U/. {k — 1, 2, 3) and we therefore have, 
neglecting the physically unimportant constant term: 

7 3 , v i 

t^r.M,*) (-) ( 12 ) 

1=1 fc=i \ z o/ 

where zq is an arbitrary length, chosen to be characteris- 
tic of the trap geometry, and set equal to 1 mm for trap 
A. The parameters Vki only depend on the geometry of 
the electrodes and are determined from our fit on the nu- 
merically calculated grid. Due to symmetries one can set 
Vxi — for I > 2, V21 — for I odd and V31 — for I 
even. The results of the fits are presented in Tabled The 
relative difference between the results of SIMION and the 
fit is smaller than 1% within a radius of 400 /im from O. 

We have also designed another trap geometry (trap 
B) which achieves tighter confinement and smaller atom- 
clcctrodc distances [see Fig. 0c)]. Such a setup is neces- 
sary for atom-surface interaction studies. Trap B is de- 
rived from trap A by bringing the ring of the top plane 
of Fig. 0b) down onto the lower plane. The electrode size 
is reduced by a factor of 10 (and zq accordingly reduced 
to 100 /im) while the distance between the two planes 
remains unchanged at 1 mm. Numerical calculations with 
SIMION show that there exists a saddle point for the elec- 
tric field modulus at a point O, 120 /im above the lower 
plane surface. It is possible to apply the same numerical 
treatment to this trap as before and fit the electric po- 
tential around O (here however we have fewer symmetries 
and can consequently set fewer of the vm equal to zero). 
This procedure produces a fit correct to 0.1% within a 
radius of 50 /im of O. The results are shown in Tabled 



2.4 Simulation of trajectories 

The evolution of the atom in the time- varying electric field 
being adiabatic (see Sec. 12. 1|) we can simply write £(t) — 
£[E{t)\. From V mt , it is easy to derive a formula for the 
electric field modulus E, and therefore for the potential 
energy of a trapped atom. It is then simple to derive an 
expression for the force on the atom along each direction. 

In the case of trap A, due to its symmetries, the poten- 
tial energy considered as a power series in \r\ has a small 
number of terms that dominate. The first non-trivial term 



is a quadrupolar term, exactly equivalent to Q 
3V2lUiU 3 (t)v n v 22 



£ 



Quad 



4-7T 



(2z 2 -x 2 -2/ 2 ) 



(13) 



If, once again, we have U${t) — U30 cos(wi) then for each 
direction the associated Mathieu equation will have a char- 
acteristic parameter [i = x,y, z) equal to: 



q x = q v = a- 



3V21 UJJJ30V11V33 
9 4 ' 



(14) 
(15) 



The second non-trivial term is a term linear in z, exactly 
equivalent to Ijllfl : 



(16) 



15 U1II2VHV22 



allowing us to calculate the value of U2 necessary to com- 
pensate gravity. 

The use of only these first two terms fits the exact po- 
tential energy to better than 1% within 100 /im of O and 
allows one to find solutions for the trajectories via anal- 
ysis of the Mathieu equation. This treatment, however, 
does not suffice for trap B due to its non-zero values of 
V21 and U32. Therefore, to treat trap B, and to achieve a 
more precise calculation of the trajectories, we retain a 
full numerical approach: 

— Each trajectory is computed, using an adaptive Runge- 
Kutta algorithm. 

— The total time over which we calculate the trajectories 
ranges from a few ms to several seconds according to 
the application. 

— A trajectory is considered trapped if the atom remains 
within the sphere of validity of the numerical fit of the 
potential (radius of 400 /im for trap A, 50 /im for trap 
B). 99% of untrapped atoms leave this zone within 
0.55 s. 

— We use initial conditions for the velocity and position 
of the trapped atom coherent with excitation from a 
trapped cloud of Rubidium atoms (trapping frequency 
1 kHz along each direction) at thermal equilibrium 
(temperature To = 300 nK) and with the loading mech- 
anism that we consider in Sec. 13.21 These correspond 
to Gaussian distributions of widths 5.35 mm/s and 
0.27 /im for the velocity and position respectively. To 
the initial velocity is added the single recoil velocity 
v r = 6 mm/s along the direction Oy that would be 
received on efficient adiabatic excitation to the Ryd- 
berg state. This temperature and trap frequency cor- 
respond to an atomic cloud close to condensation, or 
already condensed (average phonon number n p h « 6). 
The initial conditions that we have taken remain valid 
however, as long as the number of atoms in the trap is 
low enough for mean field interactions between ground- 
state atoms to be neglected. 

— For each trajectory we record periodically the atomic 
position and velocity, as well as the value of the electric 
field modulus and its direction. 
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Table 1. Results of the fit of the potential coefficients vm for trap A (zo = 1 mm). The first two lines give values for the fit 
of applied voltages Ui and Ui respectively, the third line for applied voltage U3 with rj = 1. It is possible to cancel the electric 
field created by C/3 at the centre of the trap by setting 77 = 4.49 (fourth line) 
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Table 2. 

Re- 
sults 
of 
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fit 
of 
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po- 
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co- 
ef- 
fi- 

cients 

vu 
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(20 
()fj /j,m). 
The 
dif- 
fer- 
ent 
lines 
have 
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mean- 
ing 
as 
in 
Ta- 
ble!]] 



In all cases we check in the shallow electrodynamic trap 
that the trajectory extension, of a much greater size than 
the initial, tightly confined ground-state cloud, is also much 
larger that the typical de Broglie wavelength of the atom 
(of the order of 0.50 /im). This justifies a classical calcu- 
lation of the motion. 



2.5 Trap performance 

We have studied the trapping characteristics for a wide 
range of parameters, especially the voltage U\ which con- 
trols the mean electric field experienced by the atoms, 
the voltage C/30 and the frequency ui of the a.c. field, the 
last two together controlling the confinement. For all of 
these parameters we find qualitatively the results pre- 
dicted by the Mathieu equation. Bound trajectories are 
found to be composed of two motions: a fast micromo- 
tion, of frequency uj and of relatively small amplitude, 
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Fig. 3. (a) A typical trajectory of an atom in trap A, shown 
in the xOz plane. The atom is initially at the origin of the trap 
and has a velocity (12,0,12) mm/s. (b) x (solid line) and z 
(dotted line) coordinates for the same trajectory as a function 
of time. 



and a slower macromotion of larger amplitude and of fre- 
quency uj x = ujyrti u> z /2, an order of magnitude slower 
than lj (see Fig.yJ). In the case of non-bound trajectories, 
the amplitude of the macromotion usually increases slowly 
until the atom is finally expelled from the region in which 
Vf it is valid. 

The trapping efficiency as a function of the frequency 
u> is obtained by simulating 100 trajectories and record- 
ing the fraction of them remaining trapped after 1 s. The 
results for trap A are shown in Fig. 0] No trapping is 
possible below a threshold frequency uth/^ir. This is in 
qualitative agreement with the prediction of the Mathieu 
equation stability criteria \q z \ < 0.907. For the parameters 
of Fig. 0| this would set a threshold frequency at 395 Hz, 
compared to the value of cj^/27r =400 Hz found in the 
simulations. Just above u>th the efficiency depends essen- 
tially on the initial temperature To of the sample and can 
be very close to 100% for comparatively cold atoms. If 
one continues to increase u> one observes a slow decay of 
the efficiency. This can be easily understood as due to the 
increasing 'averaging' of the time-dependent electric field 



| 800 1600 2400 3200 4000 4800 

CO, h /27t 

OV27I (Hz) 

Fig. 4. Trapping efficiency (trap A) as a function of frequency 
uj. Ui =0.2 V, U 2 = -3 mV, U 30 =0.056 V, T =l yK. 



seen by the atom. The Mathieu equation predicts that the 
macromotion frequencies u> x , lu v and w z decrease as 1/lo. 
There will therefore rapidly come a point when uj is so 
much larger than u) x , oj y and lo z that the time-dependent 
field Es(r,t) averages completely and the atom, seeing a 
homogeneous field, is no longer trapped. 

In the following, uj will be chosen slightly above LOth 
where the trapping efficiency is largest. Although this does 
not correspond to a situation where a clear cut separa- 
tion of the micro- and macromotions exists, we can still 
distinguish between the two. We present in Table |21 the 
frequencies of the macromotion oj p = lo x = oj y and ui z for 
different geometries and sets of voltages. In order to es- 
timate the depth of each trap we calculate the trapping 
efficiency for different initial temperatures. We define Td 
as the temperature above which more than 50% of the 
atoms are lost after 1 s. For trap B, two values of U\ are 
shown (0.2 and 1.5 V), corresponding to two different bias 
electric field amplitudes. Comparison of the results clearly 
show that the smaller the trap and the larger the bias 
field the tighter the confinement. Trap depth can reach as 
high as the mK with motion frequencies in the kHz range. 
However, most of the following simulations have been per- 
formed with a small voltage C7i, necessary for studies of 
coherent superpositions of \g) and |e), as we shall see in 
Sec. [5] Even in these unfavourable conditions for trapping, 
the depths attained are compatible with typical cold-atom 
temperatures and oscillation frequencies remain reason- 
able, around 100 Hz. 

Small departures of the electric field direction from 
u z in the trapping region, due to transverse components 
of E2 and E3, will later prove critical in calculating the 
efficiency of the spontaneous emission inhibition (Sec. 0} 
and the rate of dephasing between \g) and |e) (Sec.|SJ). We 
therefore note in Table [3] the mean angle (9) between the 
local electric field and u z experienced by an atom over the 
course of its trajectory. 
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Type 


V 


L7i (V) 


U 30 (V) 


U 2 (mV) 


w/2tt (Hz) 


uj p /2-k (Hz) 


u! z /2k (Hz) 


T d (/iK) 


(0) (mRad) 


Trap A 


4.49 


0.2 


0.056 


-3 


430 


64 


175 


180 


9.38 


Trap B-a 


0.05 


1.5 


0.5 





20700 


1460 


2910 


1000 


4.22 


Trap B-f3 


0.05 


0.2 


0.14 


-0.45 


2860 


207 


414 


35 


3.57 



Table 3. Trap A and trap B correspond to the geometries of Fig.|2b) and (c) respectively. Labels a and f3 for Trap B correspond 
to different voltage settings. For each trap, we give the values for U\, U2, U30 and the frequency lu/2-k. The frequencies uj z /2tt 
and Up/2-K are the typical longitudinal (along Oz) and transverse oscillation frequencies in the trap. The trap depth Td is defined 
as the temperature for which half of the atoms remain trapped. (8) is the mean angle between the local electric field E and u z , 
averaged over 100 trajectories at temperature T d /2. 
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Fig. 5. Level diagram demonstrating the dipole-blockade ef- 
fect. 



3 Loading of the trap, preparation of a 
subpoissonian sample of Rydberg atoms 

3.1 The dipole-blockade efFect 

For many quantum information processing experiments, 
the deterministic preparation and storage of an individ- 
ual qubit is a critical requirement. This has already been 
achieved, for example, with trapped ions. In previous ex- 
periments with circular Rydberg atoms [HJ, we expressly 
set the rate of excitation extremely low in order to be 
able to neglect two-atom events. This is, of course, at the 
expense of very long data acquisition times. We propose 
here to deterministically load the trap with a single Ryd- 
berg atom by using the dipole-blockade effect. This phe- 
nomenon has been proposed as an efficient way to perform 
quantum gates between two Rydberg atoms |13j but can 
also ensure the preparation of one and only one low an- 
gular momentum Rydberg state by laser excitation of a 
dense micron-sized cloud of ground-state atoms [201 • The 
essential idea behind this phenomenon is that level shifts 
due to the dipole-dipole interaction make the laser nonres- 
onant for the excitation of more than one Rydberg atom. 

In order to assess the effect of atom- atom interactions 
we consider an ensemble of N cold ground-state atoms 
(temperature in the [iK range) in the presence of a laser 
resonant with the transition between a low- lying state \b) 



and a given Rydberg state |r). Atomic motion can be ne- 
glected during the laser excitation and hence each atom 
position Ri is treated as fixed. Let us assume that the ini- 
tial trap confines the atoms in a micron-sized region such 
that \Rij\ = \Rj — Rj\ < 1 /Um for any couple of 
atoms. We then switch on the laser light, of frequency ojl. 

Let us consider for a moment only atoms 1 and 2. 
Fig. shows the relevant levels for this simplified situa- 
tion. We tune the laser frequency to be resonant on the 
transition from the lower state \b1b2) to the singly excited 
state, where one of the atoms has been excited to the Ry- 
dberg state |r). In the absence of interactions, this laser 
would also be resonant with the transition towards the 
state |rir2) where both atoms 1 and 2 are excited into the 
state |r). We must, however, take into account the cou- 
pling of this level \rir 2 ) with other doubly excited states 
of different Rydberg levels, for example Ir^'), shifted in 
energy from \r1r2) by 5dd (see Fig. [5J. This coupling, of 
magnitude Wdd, will cause a shift Add of the level |rir 2 ) 
which, if large enough, will make the excitation of the 
doubly excited level nonresonant and therefore strongly 
suppressed. We must therefore calculate Add- First order 
perturbation theory gives us: 



A wj d 

^dd = -J • 

Odd 



(17) 



Sdd will be of the order of the separation of Rydberg levels, 
a few GHz for n ~ 50. A classical order of magnitude for 
W dd is given by: 



hW A . 



87re |i2 12 | 



(18) 



where ao is the Bohr radius. For a distance I-R12I = 1 M m 
and n = 50 one calculates W dd = 3 GHz. We therefore 
calculate a shift of the doubly excited level Add ~ 1 GHz. 
This order of magnitude calculation is in agreement with 
the result Add = 100 MHz found by Saffmann and Wal- 
ker [20| • This shift is easily resolved using modern lasers 
and hence excitation of multiple Rydberg atoms will be 
out of resonance with the excitation laser. It is shown 
in [201 that for the resonant excitation of the transition 
1 6) — > |r) in a cloud of ground state atoms of diametre 
5 /jm using a laser of Rabi frequency 1 MHz one obtains 
multiple- or non-excitation of |r) for less than 1 event in 
10 (for a number of ground-state atoms of the order of 
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a few hundred, compatible with the assumption that we 
can neglect interactions in the calculation of the initial 
conditions). 



3.2 Loading from a magnetic atom-chip trap 



Atom-chip traps can fulfill all the conditions required 
by our proposal. They allow for very high densities of 
ground-state atoms in trap volumes as small as a few 
fim. 3 |3T]. They can, by construction, store atoms close 
to surfaces and the operation of a conveyor belt, ideal for 
bringing atoms from a capture region into the Rydberg 
trap, has been proven [221 - Finally, evaporative cooling in 
such traps allows one to attain temperatures as low as 
a few hundred nK [231 - The need for operation at cryo- 
genic temperatures, in order to get rid of blackbody radi- 
ation which would rapidly destabilize the Rydberg states 
(see Scc.0J, makes it impossible however to use standard 
atom-chip designs with normal conductors. We plan to 
use instead superconductors to create the magnetic fields. 
This would allow us to pass current without Joule heating 
as long as one remains below the critical current. In stan- 
dard atom-chip experiments at room temperature the cur- 
rent densities are typically of the order of 5 x 10 6 Acm~ 2 
123] • We have checked that it is possible to reach similar 
current densities with our superconducting wires. In the 
case of a superconducting Niobium wire of thickness 1 /itm 
and width 10 /im, sputtered on a Silicon Oxide substrate, 
we find I c — 0.5 A, corresponding to a current density of 
5 x 10 6 Acm~ 2 . 

Thanks to our microfabricated design, it is simple to 
integrate our Rydberg trap with the magnetic atom-chip 
trap and also the electrodes necessary for the transfer of 
the atom from a low angular momentum Rydberg state 
to a circular state [2HJ - We have designed a multi-layer 
configuration with the electrodes necessary for the Ryd- 
berg trap placed on top of the wires of the magnetic atom 
chip. The circularization is achieved using wires on the 
magnetic atom-chip layer to create the necessary mag- 
netic field along with the application of a r.f. voltage to 
electrodes on the Rydberg trap layer. The Rydberg trap 
electrodes will not significantly affect the magnetic field 
created by the wires of the lower layer as long as a normal 
conductor such as Gold is used. Evaporative cooling in 
the magnetic atom-chip trap can provide a few hundred 
atoms in a micron-sized cloud at temperatures as low as 
a few hundred nK [23], conditions ideal for the operation 
of the dipole-blockade effect. The magnetic trap is then 
suddenly switched off and the excitation towards a low an- 
gular momentum Rydberg state performed immediately. 
Carefully designed laser excitation schemes result in the 
atom receiving only one optical photon recoil, just before 
the circularization process, itself lasting about 20 fxs [23 
[B] . The flexibility of magnetic atom chips will allow us to 
superimpose the centre of the magnetic trap with that of 
the Rydberg trap. 
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Fig. 6. Ideal geometry for spontaneous emission inhibition. 

4 Making the atom long-lived by spontaneous 
emission inhibition 

4.1 Principles 

The basic idea behind the inhibition of spontaneous emis- 
sion is to place the atom inside a cavity containing no 
mode into which the atom can emit its photon upon de- 
cay [21] • For the sake of simplicity we first of all consider 
the geometry of Fig. [5] where a circular Rydberg atom is 
placed between two infinite, perfectly conducting planes 
separated by a distance L. This situation is relatively close 
to that of Fig. [21(b). We now assume that there exists be- 
tween the two planes a d.c. electric field E dir orthogonal to 
the planes (and hence parallel to their normal u z ). Obey- 
ing the standard selection rules, a circular Rydberg state 
of principal quantum number n possesses only a single 
possible decay channel, this being towards the lower cir- 
cular state n— 1 (for example |e) — > \g)). In the process it 
emits a photon of polarization a + with respect to u z , of 
wavelength X sp (X sp ~ 6 mm for the decay of \g) and |e)). 
Due to its polarization, the electric field associated with 
this photon is perpendicular to u z and hence parallel to 
the planes. The maximum wavelength of any such mode 
inside this cavity is 2L as the electric field must cancel at 
the surface of each plane. Therefore, if A sp > 2L, there 
exists no cavity mode at resonance with the spontaneous 
emission transition and the atom cannot decay, remaining 
in the circular state for an infinitely long time. 

We now detail a more rigorous, semi-classical approach 
to this phenomenon, as found in Ref. |2U- We calculate 
the effect of the field radiated by the dipole upon itself, 
and hence the consequent energy shifts and lifetimes. The 
electric field radiated by a dipole D — Dexp(— iu>t)u d at 
position r d oscillating at frequency u> along unit vector u d 
is given by: 



E di p(r,t) = 



Duo 6 



(19) 



where ^ cav "> (r, r d , u>) is the susceptibility of the field in- 
side the cavity, characterizing the linear response of the 
system to the dipole movement. This susceptibility can 
be written in the form: 

7*-™* (r, r d ,w) = F i0) (r, r d , w) + T {r) (r, r d , w). (20) 

The field corresponding to ^°\r, r d , uS) is the field radi- 
ated by the dipole in free space and that corresponding to 
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T^ r ' (r, r<j, to) reflects the modification introduced by the 
presence of the cavity and is easily understood as being 
the field radiated by the images of the dipole in the sur- 
faces of the cavity. When we calculate the interaction of 
the dipole with the total field we will therefore have two 
terms: 

— A first representing the interaction of the dipole with 
its own free space field. This term leads us to the cal- 
culation of the natural lifetime of the dipole and the 
Lamb shift of its frequency. This calculation however is 
plagued by divergences. Their values must be derived 
by a fully quantum treatment, and we simply insert 
these results into our calculation. 

— A second corresponding to the interaction of the dipole 
with the field of its images. This leads to the modifi- 
cation of the natural lifetime (that we search to calcu- 
late) along with the accompanying frequency shift. Its 
calculation, in contrast to the first, poses no problem. 

This approach gives the following expression for the 
rate of decay of the dipole (atom) inside the cavity |2fi| : 



r = r 1 



-Im 



(21) 



the imaginary part of the susceptibility corresponding to 
the component of the image field in quadrature with the 
dipole oscillation and hence responsible for dissipation. We 
present in Fig.[7Jthe evolution of the ratio r/To as a func- 
tion of the distance L in two specific situations: D parallel 
to u z ; D perpendicular to u z . In the latter, the case of a 
emission in the presence of a quantization axis along u Zl 
one observes a complete inhibition for L < A/2. In the 
former, the case of ir emission in the presence of the same 
quantization axis along u Zl the dipole couples to modes 
obeying different boundary conditions and its emission is 
exalted for small L values. In the case of a a transition 
in the presence of a quantization axis slightly tilted with 
respect to u z we have components of D both parallel and 
perpendicular to u z and consequently the possibility of 
decay via both inhibited and exalted channels. Inhibition 
has already been observed in the microwave domain [27J 
as well as for optical transitions [2H| • A value of less than 
0.1 for r/To has been reported in the case of Rydberg 
atoms. This would already increase the radiative lifetime 
of circular states to a significant fraction of a second. 



4.2 Limitations of the inhibition 

The ideal case presented in Sec. 14. H is, however, far from 
the realistic geometries of Fig. Etb) and (c), the imper- 
fections of which we must now evaluate. As a first exam- 
ple, the electrodes do not form perfect infinite mirrors. It 
is possible to neglect the effect of their finite size if, as 
will be the case, the transverse dimensions of the trap are 
much larger than X sp . A more critical effect comes from the 
residual absorption by the imperfect mirrors, dissipating 
the energy radiated by the atom and limiting its lifetime. 
This effect must therefore be included in the calculation 



r 2 




Fig. 7. Spontaneous emission rate F, in units of Fo, as a func- 
tion of the cavity plane separation L. The dipole is either par- 
allel to the mirrors (r« , continuous line) or perpendicular to 
them (F±, dashed line). Both are calculated for the case of 
perfectly conducting mirrors. 



of the susceptibility T (r '{r, r d , u). The dipole of an image 
produced by a number of reflections in the mirrors Nr is 
therefore multiplied by a factor: 



(pe«) 



(22) 



where \p\ 2 is the reflection coefficient of the mirror and x 
accounts for the phase shift at each reflection. For trapping 
electrodes made of ordinary conductors (such as Gold) p 
and x are related to the skin depth 5 by the following 
expressions: 



P 

X 



2ttS 
ex P( 

2tt5 



(23) 
(24) 



provided that, as is the case here, 5 <C A • Figure^a) 
presents the inhibition factor -Til /.To as a function of the 
skin depth 5 in the case of a dipole radiating at 50 GHz 
(A = 6 mm). The cavity spacing L is 1 mm and the dipole 
120 /xm away from one mirror as in trap B. One can see 
that inhibition is still efficient (a factor of greater than 
100) for skin depths of less than 100 nm. Operating, as 
we must, under cryogenic conditions, we can therefore use 
Gold, for which (5 = 30 nm at 1 K. For imperfect mirrors 
there is another effect that must be considered. The inhi- 
bition factor now depends upon the position of the atom 
between the two planes. Figurc|H{b) shows this effect. The 
closer to the surfaces, the worse the inhibition. One might 
fear that this would pose problems in trap B. However, 
even in the case of a circular Rydberg atom in trap B, 
120 fj,m away from Gold trapping electrodes at IK, the 
inhibition factor of 340 would imply a radiative lifetime of 
10 s. 

We must also consider the consequences of our quanti- 
zation axis, following adiabatically the electric field E(r, t), 
being not always perfectly perpendicular to the plane mir- 
rors. A non-zero angle 9(r, t) between E(r, t) and u z will 
be responsible for a small component D± = D sin 9 of the 



John Mozley et al.: Trapping and coherent manipulation of a Rydberg atom on a microfabricated device: a proposal 11 



(a) 

8x10-3 
o 6x10-3 
^ 4x10-3 
2x10-3 




20 40 60 80 100 




0.1 0.2 0.3 0.4 0.5 

z (mm) 



whatever the distance between the dipole and the mir- 
rors. Under the conditions considered (Table OJ), r colI /Tq 
remains greater than 150. 

The final effect that must be considered in evaluating 
the residual lifetime of a circular Rydberg atom in the trap 
is excitation by blackbody photons in the modes of the 
cavity. These can induce 7r-polarized transitions between 
\g) (resp. |e)) and the other states of the m — 49 (resp. 
m = 50) multiplicity (for example \g) — > see Fig. [IJ. 
The rate of this process is proportional to the mean photon 
number n t at the relevant transition frequency times the 
relaxation rate r v for each transition, the total rate -T b .b. 
being equal to the sum of all possible transitions -Tb.b. (k) = 
Y,k> n t {kk')r^{kk'). We calculate r B . B .(g) = 3.15 s^ 1 and 
r B . B .( e ) = 2.75 s -1 , the effect of exaltation of it transitions 
having been taken into account. If the trap is cooled down 
to 1 K, one has nt = 0.1 and a corresponding radiative 
lifetime of about 3 s, once again not a limiting factor. 



5 Preserving atomic coherences over long 
times 



Fig. 8. (a) Spontaneous emission inhibition factor Fm/Fo as 
a function of the skin depth 8. The conditions are the same 
as in Fig. |S| with L — 1 mm, A = 6 mm and a dipole 120 /im 
away from one mirror, (b) Fn/Fo as a function of the atom- 
mirror separation for the same values of L and A as (a) and 
with <5 = 30 nm, corresponding to the case of Gold at 1 K 



dipole being orthogonal to the mirrors (and for a reduced 
dipole D|| — DcosO parallel to the mirrors). It is there- 
fore possible for the atom to couple to modes possessing 
exalted decay. Taking into account only this decay chan- 
nel, the probability for remaining in the excited state is 
therefore given by: 



P = cxp f- jT r ± sin 2 9 [r(t),t] dt 

= exp (~r J j± sin 2 9 [r(t),t] dtj , (25) 



where the integral is performed along the trajectory fol- 
lowed by the atom. 

As we have seen (Tabic |3J), *C 1 for the trajectories 
considered and we can therefore approximate Du = D 
and sin 2 9 = 9 2 . From Eq. EH one can deduce a corrected 
inhibition factor in r corx /Tq, related to the value of 9 along 
a trajectory: 

% = 5l + ^ 



(26) 



where: 



2 = — / f)-,ll. 

ttr. 



(27) 



is the average of the square of the angle 9 over the total 
duration of the trajectory, tt r .- For L = 1 mm, A = 6 mm 
and (5 = 30 nm we find an exaltation factor r±/To of 4.5, 



To this point, we have shown that it is possible to trap a 
single atom in a circular Rydberg state over times in the 
second range. It is, however, important to see if one can 
manipulate and maintain its internal state with the same 
precision. More precisely, if we hope to use this technique 
in quantum information experiments or for high precision 
spectroscopy, it is important to verify that one can prepare 
an atom in a coherent superposition of states and probe 
its phase at a later time. In the following we will focus on 
the transition between the two circular states \e) and \g). 



5.1 Electric shifts as the main cause of dephasing 

As we saw in Sec. 12.11 the Stark polarizabilities of levels 
|e) and \g) are slightly different. Fig. [5] shows the energy 
levels as a function of the electric field amplitude. The 
frequency v eg = {£ e — £ g )/h of the transition has a strong 
dependence on the electric field: 



5v eg = v eg {E) - v eg (0) = 5a eg E 2 



(28) 



with 5a eg = -25.5 Hz/(V/m) z . If one considers the case 
of an atom in trap A in the conditions of Table prepared 
from a cloud at 0.3 /iK, it experiences over its trajectory 
a mean electric field amplitude E a = 400 V/m, with an 
excursion of AE — ±1 V/m. The corresponding frequency 
broadening: 

Av eg = 2Sa eg E a AE, (29) 

is of the order of 20 kHz. This broadening is inhomoge- 
neous as it is different from one trajectory to the next and 
could only be controlled by perfect control of the initial 
conditions of the atomic motion. The dephasing time as- 
sociated to this motional broadening is of the order of a 
few tens of us. 
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In addition to this dephasing, the trapping frequencies 
(u p , uj z ) for states |e) and \g) differ by about 10%. The 
trajectories for the two states are therefore rapidly sep- 
arated ('Stern Gerlach' effect). Coherence would be lost 
when this separation exceeds the wave-packet coherence 
length, of the order of the de Broglie wavelength (about 
0.5 /Ltm for the conditions considered). 



5.2 Tailoring the atomic levels to cancel electric shifts 

A similar situation is observed for ground-state atoms in 
magnetic or optical dipole traps. The potential of these 
traps is also, in general, level dependent. However, well 
chosen trap laser wavelengths [3U| or bias magnetic fields |31l 
I32j minimize the effect. To achieve the same end in the 
Rydberg trap we propose to use a microwave field dress- 
ing in order to tailor the atomic energies. A 7r-polarized 
microwave field, of frequency uo, is fed inside the trap. 
This field is compatible with the trap boundary condi- 
tions. For the sake of clarity, and as alluded to earlier, 
we consider that additional boundary conditions, far from 
the trap center, support a propagating wave in a guided 
mode: its amplitude, independent of z and x, is maximal 
at the origin and varies sinusoidally with y, canceling at 
y = ±1 cm. 

One can easily understand the effect of the microwave 
by the consideration of a simple two-level approximation 
where the microwave only couples levels \g) and \i) (see 
Fig. HJ). The state \i) experiences a large, linear Stark 
effect [1 MHz/(V/m), see Fig. EJ. We define L? as the 
classical Rabi frequency on this particular transition, de- 
pendent on the applied microwave power. The field is 
detuned by So — uj g i(E a ) — u>q > to the red of the 
\g)/\i) transition, Stark shifted by the average field E a . 
As a result the 'dressed' level \g) acquires a fraction of 
the Stark polarizability of \i) and the dependence of its 
energy on electric field can be made parallel to that of 
|e), at least in the neighbourhood of E a . More precisely, 
we can choose the values of the two independent param- 
eters Sq and £2q so as to cancel the linear and quadratic 
terms in the expansion of the dressed state transition fre- 
quency uj e g(E) around E a . For E a = 400 V/m, the result- 
ing d = 2tt x 746.158 MHz and f2 = 2ir x 228.442 MHz 
have reasonable values. The remaining higher order terms 
in the expansion of the transition frequency uj eg (E) are 
smaller than 0.05 Hz over the complete field amplitude 
range of AE = ±1 V/m around E a . This would result in 
dephasing times longer than 20 s. 

We note here that this scheme would also compen- 
sate for the problem of patch effects on the surface of the 
electrodes if the field which they create remains relatively 
small compared to the directing field amplitude E a . This 
effect is caused by the misalignment of crystal orientation 
in adjacent crystals in the metal. Let us consider a sur- 
face with voltage difference AV between adjacent crystals 
of characteristic size a. An ensemble of such patches will 
create at a distance d above the surface a field: 
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Fig. 9. Energies of levels \g), |e) and \i) as a function of the 
electric field. States \g) and e) have a small, approximately 
quadratic Stark effect, magnified in this figure by a factor 5 
for the sake of clarity. Level \i) has a much larger, linear Stark 
effect. The microwave dressing (vertical arrow) predominantly 
mixes \g) and \i) and reduces the difference in Stark polariz- 
ability between |e) and \g), the resulting dressed level. 
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Patch-effects for Gold electrodes were measured at much 
closer distances in Ref. |33) where they found E pEltch ~ 
10 4 V/m for d = 0.25 /im and a ~ 30 nm. Being conserva- 
tive, we assume the patch size to be 100 nm. Building on 
these figures we estimate the field to be 0.008±0.0003 V/m 
for trap A and 0.15 ± 0.005 V/m for trap B, the inho- 
mogeneity being calculated over a volume of dimension 
Ar — 15 /im and 5 fim, the typical excursion of atoms 
in traps A and B respectively. All of these values remain 
much smaller than E a and would hence be well compen- 
sated for by our technique. 



5.3 Assessment of performances 

The simple explanation of Sec. 15.21 is. however, far from 
describing the true physical reality: 

— The value calculated for So is not significantly smaller 
than the energy difference between state \i) and state 
(see Fig. 0, which is therefore also significantly 
coupled to \g). Moreover the microwave field also cou- 
ples state |e) to the n = 52 manifold. It is therefore 
important to adopt a multilevel approach in order to 
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account for the coupling by the 7r-polarized dressing 
microwave of all the levels of equal to (see Fig. 

— This coupling therefore mixes the level \g) with other 
levels of the same to — 49 (as is its aim) and equally 
|e) with levels of m = 50. It must not be forgotten, 
however, that these levels have finite, potentially short, 
lifetimes in the cavity. We can therefore now see a form 
of spontaneous emission by \g) and \e) in which they 
absorb a photon from the dressing microwave mode 
and immediately re-emit it into another mode of the 
cavity. These spontaneous events would destroy any 
coherence existing between \g) and |e). We must there- 
fore quantify the rate at which these events occur and 
check that the corresponding lifetime is not limiting. 

— The dressed level transition frequency uj eg is dramati- 
cally dependent on the microwave Rabi frequency, pre- 
viously considered constant and equal to J?o- But Q 
is actually time- and position-dependent due to the 
spatial profile of the mode (nodes at y = ±1 cm, see 
Sec. 15.2(1 and due to the variation of 9(r,t). The in- 
stantaneous Rabi frequency experienced by an atom 
is therefore given by n(r,t) = fl f{r) cos 6(r 7 t). This 
variation of the Rabi frequency creates an inhomoge- 
neous broadening which must be included in our final 
calculation. 

— The tilt also couples a fraction of the microwave 
power to a transitions within the Rydberg levels, with 
an effective Rabi frequency n a (r, t) cx S7of(r) sm9(r, t), 
significantly smaller than fl. The large detuning Sq pre- 
vents us from having resonant one-photon transitions 
between the levels |e) and \g) due to this coupling. 
Nevertheless, there could exist a multiphoton transi- 
tion coupling the circular states to adjacent manifolds. 
A 1 Hz Rabi frequency would be sufficient to transfer 
population from the circular state to this adjacent one, 
and limit the coherence time between |e) and \g) to 1 s. 
We must therefore check that no such resonance exists. 

— To achieve the necessary precision in the position of 
the levels |e), \g), \i) and \i') we must go beyond the 
quadratic Stark expansion of Eq. ^ 

We will now address these points one by one, turning 
first of all to the last. To achieve the precision necessary 
we diagonalize the Stark Hamiltonian of the circular state 
and the 4 manifolds of greater n above it, in the presence of 
an electric field. The resulting energies were fitted by a 4th 
order polynomial over electric fields from to 1000 V/m 
with an error less than 0.5 Hz. It was checked that adding 
a 5th manifold above the circular state to the Hamiltonian 
changed the resulting energies by less than 1 Hz. 

In order to take account of the coupling of levels other 
than \g) and \i) by the microwave field, we must carry out 
a diagonalization of a Hamiltonian containing all relevant 
couplings. The microwave field is fully defined by its classi- 
cal Rabi frequency at trap centre on the \g) — > \i) tran- 
sition and its detuning <5o- The angle 9(r,t) being small, we 
consider firstly only the 7r-polarized component of the mi- 
crowave field, coupling levels of equal to. We have adopted 
a dressed level approach to this problem and therefore 
consider the ladder of levels InjU^m,^) where the ad- 



ditional quantum number counts the number of pho- 
tons in the microwave field. For the case to = 49, relevant 
for level |<?), the coupling between initial and final states 
\rii, ni^i, 49, + 1) and |n/, rt\ j, 49, n T ) is given by: 

K) = ?Znli V^+T (32) 

where T n f t n±f is the Rabi frequency on the considered 
transition for — 0, proportional to the dipole matrix 
element A4 n {'nt'i ■ However, as the dressing microwave is a 
classical field (and hence n v 3> 1), one can assume y/n^ « 
y / n^~+T « constant. Therefore each coupling strength can 
be normalized to according to: 

JAu.ni.i \ n K) - U 51|1 - J/ ■ .51,1 • ( 66 ) 

We were therefore able to construct the matrix of a Hamil- 
tonian at constant m including M multiplicities in n (ie. 
n = to+1, • • • , to + M), and 2A+1 multiplicities in (ie. 
n,,- — N, • • • , n v , • • • , n-n + N) . For the initial level positions 
£(k;E) we use the quadratic Stark expansion of Eq. ^ 
for all levels apart from |<?), |e), \i) and \i'}, as explained 
above. After diagonalization of this Hamiltonian, we are 
left with a new ladder of levels, £(k;E). These levels are 
slightly shifted in energy with respect to the unperturbed 
states and it is of course this shift that we hope to use in 
order to render £{g,n^\ E) and £(e,n 7r ;£') parallel. This 
shift being small, and there being no resonant coupling 
due to the dressing field, we are able to associate each of 
the dressed levels unambiguously with one of the unper- 
turbed levels. We label the state arising from \n, rii, m, n n ) 

as \n, tlx, to, Utt) where the new quantum numbers charac- 
terize the now mixed atom-field state. We computed the 
change in energy of £{g, E) and £(e, n^; E) as one in- 
creases M or N by one, hence increasing the size of the 
Hamiltonian diagonalized. We found this difference to de- 
crease exponentially and fall below the Hz level for M > 6 
and N > 4. By a diagonalization of the Hamiltonian with 
M = 6 and N = 4 we were therefore able to calculate 
£(g,n„; E) and £{e,n^] E) to the Hz level for a given 
electric field E, microwave power (corresponding Rabi fre- 
quency i?o) and detuning 5q. We can therefore write: 

£ k (E, Q ,So) = £(k, E) - n^hujo, (34) 

for the potential energy of a trapped atom in state k = g 
or e, in the presence of the microwave dressing field. 

From these energies one can determine the dressed 
level transition frequency: 

Ueg{E, n 0) S ) = [£ e (E, Ho, 5 ) - £ g (E, fl , S )] /h, (35) 

and analyse its dependence on the electric field. In the 
neighbourhood of E a we have: 

CJeg(E, n ,5 Q ) = UJeg{E a , ^0, S ) (36) 

+L(E a ,n ,S Q )(E-E a ) 

+Q(E a ,n ,S a )(E-E a ) 2 

+0((E-E a f), 
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and it is the coefficients L(E a , J?o, So) and Q(E a , Slo, So) 
that we hope to cancel, as we did under the simpler, two- 
level treatment laid out in Sec. 15.21 Unfortunately, due to 
the effect of the microwave field upon level |e) as well as 
\g), we find that, while it is possible for any So to find 
an Qq such that L {E a , Qq, So) = 0, it remains nonethe- 
less impossible to simultaneously cancel Q(E a , (lo,So)- We 
therefore content ourselves with canceling L(E a , J?0; <^o) 
and bringing Q(E a , f?o,So) close to its minimum for the 
reasonable parameter values Qq — 2tt x 200.000 MHz and 
So = 2tt x 555.907 MHz. The remaining frequency disper- 
sion corresponding to the field excursion AE = 1 V/m is 
of the order of 10 Hz. 

Turning now to the second item on our list of imper- 
fections in the simplistic treatment of Sec. 15.21 we must 
analyse the effect of the variation of fi{r, t) due to the 
mode profile, the angle 6(r, t) or, indeed, simple technical 
noise on the microwave amplifier. By analysing the vari- 
ation of Ueg{E a , f2o + Afi, So) with Afi we were able to 
determine that the consequent inhomogeneous broaden- 
ing is less than 10 Hz for Afljfl < 2 x 10 -7 . This sets 
a tight condition on the power stability of the microwave 
source for the dressing. Both the excursion of the atom 
in the microwave field profile and the variation of 6{r, t) 
introduce a Afi smaller than this. 

The rate of spontaneous cascades of the type \g, n^) — > 

|fc, — 1) is found by calculating the dipole D(gk) be- 
tween \g, n n ) and all the eigenstates |fc, n' n ) of lower energy 
and then summing the rates of spontaneous transitions, 
proportional to the square of the dipole. The lifetimes 
thus calculated converge rapidly with increasing M and 
N towards 11.9 s for \g, n^) and 62.0 s for \e,n^) (much 
longer because of the greater detuning of the microwave 
for transitions within the m = 50 multiplicity and the 
consequently smaller contamination of |e) by other lev- 
els). These lifetimes are much longer than those already 
imposed by other limitations as presented in Sec. 14.21 and 
are therefore not an obstacle. 

Finally, it was checked that there were no resonances 
between the dressed circular states and adjacent levels 
of different m due to the small er-polarized component 
of the microwave field. No transition was found towards 
the multiplicities of Am — ± 1 with a detuning of less 
than 100 MHz. These transitions consequently pose no 
problem. All transitions at Am = ± 2,3... have Rabi 
frequencies -C 1 Hz, these processes being of second or 
higher order, passing necessarily by a virtual transition at 
Am = ± 1 which we have just seen to be hugely nonreso- 
nant, and hence strongly suppressed. 

5.4 Simulation of the coherence control 

In order to estimate the coherence decay time T-i we simu- 
lated a Ramsey interferometry experiment in the presence 
of the dressing microwave. We simulate the trajectories of 
an atom as in Sec. 12.41 with the state dependent potential 
energy given by: 

£ k {r,t) = £ k {E(r,t),n f(r)cos6(r,t),5 ), (37) 
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Fig. 10. Phase acquired by an \e)/\g) coherence over 40 ms. 
The trajectory considered is in trap A, under the conditions of 
Tableland at T = 300 nK. 

where the £k(E, f2,S) is a series expansion of the £k to 
fourth order in E and fl around (E a , f2o,So)- E(r,t) and 
0{r,t) are the electric field amplitude and tilt, and f(r) 
accounts for the dressing microwave spatial mode (see 
Sec. EH). 

We assume that the atom undergoes two short mi- 
crowave 7r/2 pulses, resonant on the \g) /|e) transition (fre- 
quency ujeg(E a , flo, So)) at times U and i/, separated by 
a long delay T — tf — ti. For each state we follow the 
state dependent trajectory. For atoms taken from a sam- 
ple at 0.3 /xK, the average separation between the two is 
~ 20 nm, much smaller than the de Broglie wavelength 
Adb = 0.34 fim.. The Stern Gerlach effect can therefore 
be neglected. The dephasing between the two states is 
then due only to the phase (f>(T) acquired by the \e)/\g) 
coherence along the trajectory r(t): 

0(T) = [ £ e (E [r(t),t] , n f [r(t)] cos B [r(t), t] , S ) 

-£ g (E [r(t), t] , Oof [r(t)] cos 6 [r(t), t] , S ) 
-u- e ~ g {E a ,n Q ,S )]dt. (38) 

Fig- EH represents this phase evolution for a given tra- 
jectory over a 40 ms time interval. One can see that, over 
times longer than the trap oscillation period, the evolution 
is quite linear. On times scales of the order of the trap pe- 
riod, 4>{T) proceeds by small steps of amplitude A(f> which 
occur when the atom is furthest from the origin. The final 
fringe contrast for the whole set of trajectories can be eval- 
uated by C — (coscj) + sin0 ) x l 2 (where the bar denotes 
an average over all trajectories). The time evolution of C, 
shown in Fier. Illf a). provides an estimate of Ti at around 
25 ms. This is in qualitative agreement with the predic- 
tions of Sec. 15.31 The evolution of the contrast over long 
time is, however, nonexponential. After a very rapid ini- 
tial fall, C(T) decays at a much slower rate. We attribute 
this behaviour to the fraction of trajectories which remain 
very close to the center of the trap and whose phase drift 
can be much smaller than 1 over significantly long times 
[see Fig. Ullb)]. 
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Fig. 11. (a) Simulated Ramsey fringe contrast C as a func- 
tion of the time interval T between two n/2 pulses in the case 
of Trap A. Conditions are as in Table El with T = 300 nK, 
we average over 10000 trajectories. Grey line: contrast decay 
without echo. The contrast undergoes a non-exponential de- 
cay, falling to 50% in 24 ms, and to 13% in 0.5 s. Black line: 
contrast decay and then revival after a 7r-pulse at = 0.5 s, 
reaching 83% at T = 1 s. (b) The histogram of the phases ac- 
cumulated at T=0.5 s sheds light on this behaviour. A subset 
of the trajectories (white bars in the histogram) accumulates 
large phases and its contribution averages rapidly to zero. The 
other trajectories (grey bars), remaining closer to the centre of 
the trap, have a much slower phase drift and account for the 
slow decay at long times, (c) Simulated Ramsey fringe con- 
trast C(T) in Trap B under the conditions of Tableland with 
To = 300 nK, averaged over 2500 trajectories. Here we have 
applied a 7r-pulse at T = 420 ms and we see a revival of the 
contrast reaching 72% at T — 840 ms. 



As the phase acquired is almost perfectly linear with 
time for all trapped trajectories, one can try to combat 
the long term phase spreading using an echo technique 
reminiscent of photon-echo techniques. Coherence preserv- 
ing echoes have also been tested for trapped ground-state 
atoms or ions |34U35U36"] . In order to examine the efficiency 
of this technique we simulate, at a time T n after the first 
Ramsey pulse, an additional 7r-pulse on the \e)/\g) transi- 
tion, exchanging the populations of the two states. Pulse 
imperfections are taken into account for by considering, 
instead of a perfect 7r-pulse, a pulse performing the fol- 
lowing transformation: 

\g) -> cos(^)|s) + sin (^)|e), (39) 

|e) -> _sin(^±%) +cos(^)|e). (40) 

where i? can be considered as the error in the angle of rota- 
tion in the Bloch sphere, d — corresponding to the ideal 
case. The i? for each trajectory is randomly chosen from 
a gaussian distribution, centered around 0, with a disper- 
sion of 0.1 Rad. In the ideal case where $ = the phase of 
the \e)/\g) coherence is multiplied by -1 upon application 
of the 7r-pulse. During the subsequent evolution, the phase 
drift continues as before and thus returns towards zero. At 
time T = 2T n , all phases are zero to within an uncertainty 
of the order of the average phase step amplitude A<j>. 

Fig. Illf a) also presents the contrast C obtained under 
these conditions as a function of T with a 7r-pulse applied 
at = 0.5 s. We see it increase sharply up to 83% around 
T = 2T„ = 1 s. This very high contrast corresponds to an 
effective — 5.4 s. Even for atoms at To = 1 /xK, we ob- 
tain C(2Tn) =57.3%. Fig. Illf c) shows that the technique 
also works in trap B, producing an effective Ti of 2.6 s. 
More complex echo sequences can be envisaged to improve 
the final Ramsey fringe contrast. Ideal 7r-pulses repeated 
at shorter time intervals can maintain the coherence over 
time scales in the minute range, being only limited by the 
radiative lifetime. 



6 Conclusions and perspectives 

We have shown that it is possible to integrate on a mi- 
crochip the elements necessary to excite a single atom into 
a circular Rydberg state, to trap it, and to perform coher- 
ent manipulation of its internal states over seconds. The 
wealth of possible electrode geometries makes it easy to 
extend these results to the case of atomic waveguides or 
arrays of traps. We have also obtained preliminary results 
proving that it could be possible to integrate the detec- 
tion of the Rydberg states on the chip itself. This relies 
on the usual state dependent field ionization of the atom 
followed by detection of the electron by a thin supercon- 
ducting wire, such as those already used in fast photon 
detectors Our proposal would offer a complete "tool- 
box" to study and control a single quantum system in an 
extremely versatile manner. It also forms a scalable archi- 
tecture for quantum information processing. The experi- 
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mental realization of this project will require the adapta- 
tion of atom-chip techniques from room temperature to 
the superconducting regime. We are currently developing 
a cryogenic experiment able to achieve this, opening new 
perspectives for atom-chip experiments. 
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